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Abstract 

The problems of airframe structural 
dynamic representation and effects of 
coupled rotor/air frame vibration are dis- 
cussed. Several finite element computer 
programs (including NASTRAN) and methods 
for idealization and computation of air- 
frame natural modes and frequencies and 
forced response are reviewed. Methods for 
obtaining a simultaneous rotor and fuse- 
lage vibratory response, determining 
effectiveness of vibration control devices, 
and energy methods for structural optimi- 
zation are also discussed. Application of 
these methods is shown for the vibration 
prediction of the Model 347 helicopter. 

Notation 

A - airframe mobility matrix 
B - rotor impedance matrix 
El - blade bending rigidity 
r - rorce 

GJ - blade torsional rigidity 
I - identity matrix 

k - rotor frequency multiple, 1, 2, etc. 

K - stiffness matrix 
M - mass matrix 

q - airframe mode generalized coordinate 
X - airframe displacements 

- airframe mode generalized mass 

w - airframe mode natural frequency 

4> - airframe mode shape (eigen vector) 

ft - rotor frequency 

[] - matrix 

{ } - column vector 


Subscripts 


A 

- absorber, airframe 


c 

- cosine component amplitude 


H 

- hub 


k 

- rotor frequency multiple. 

1 , 2 , etc . 

n 

- airframe mode number 


o 

- zero hub motion 


R 

- rotor 


s 

- sine component 
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Superscripts 

. - velocity 
- acceleration 
t - transpose 

Prediction of helicopter airframe 
vibration involves two major problem 
areas s 

• Prediction of rotor vibratory 
hub loads 

• Prediction of airframe dynamic 
characteristics . 

The effects of vibratory hub motion on 
vibratory hub loads and effects of vibra- 
tion control devices and resulting air- 
frame fatigue stresses must also be con- 
sidered. 

Methods for independent prediction of 
vibratory hub loads and airframe dynamic 
characteristics have been developed pre- 
viously and are discussed briefly below. 
Independent determination of rotor vibra- 
tory loads and airframe vibratory response 
to these loads does not account for any 
interaction between airframe vibratory 
motion on rotor vibratory loads. One 
approximate method for accounting for 
these interactions is to assume that an 
effective rotor mass is attached to the 
airframe at the rotor hub. A more direct 
method is to compute (or measure) the 
rotor hub impedance and determine compat- 
ible vibratory hub loads and hub motions. 
This method is discussed below. A simple 
example of compatible rotor load-hub 
motion is given for a single rotor heli- 
copter with vertical hub motion. In 
addition, flight test results for the 
Model 347 helicopter are compared with 
vibration predictions obtained using a 
coupled rotor/airframe vibration computer 
program. 

Rotor Vibratory Hub Loads 

Methods and digital computer programs 
have been developed for prediction of 
rotor vibratory hub loads for constant 
speed level flight conditions i/2,3. 

Rotor blades are represented by lumped 
parameter analytical models as indicated 
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in Figure 1. Iteration techniques are 
used to compute individual blade deflec- 
tions and aerodynamic and inertia load 
distributions at integer multiples of the 
rotor rotating frequency. The total 
rotating and fixed system rotor vibratory 
hub loads are obtained by summing indivi- 
dual blade root shears and moments . The 
vibratory hub loads may be computed assum- 
ing no hub motion. If the vibratory hub 
motions are known, effects of these 
motions may be included when computing 
blade aerodynamic and inertia loads. 

Airframe Dynamics 

Structural Model, Natural Modes and Fre- 
quencies, ana Forced Response 

Finite element methods have been used 
in the helicopter industry for some time 
for prediction of airframe dynamic charac- 
teristics^. As indicated in Figure 2, 
developing a finite element airframe model 
consists of: 

• Defining nodal data 

• Defining elastic properties of 
members connecting nodes 

• Defining mass properties asso- 
ciated with each node. 


The solution for natural modes and 
frequencies is made using the reduced 
stiffness matrix and the mass matrix 
associated with the retained degrees of 
freedom. 

The airframe motions are expressed 
in terms of natural modes: 

W = M M < 4) 

and, after assuming sinusoidal motion 
with no external forces. Equation (3) 
becomes: 

(5) 


The modal generalized mass is then 
computed. A value of modal damping is 
assumed for each mode, and these modal 
properties are used to compute airframe 
response to vibratory hub loads: 

= {*n} T [ M ] {*n} (6) 

<in + 2 «n“n4n + “n^n = { < J>Rn} T { F R}/»tn < 7 > 


Nodal data and properties of struc- 
tural members are used to develop stiff- 
ness matrices for individual members. 

These matrices relate forces at each node 
to nodal displacements. The stiffness 
matrices for individual members are super- 
imposed to obtain the stiffness matrix 
for the entire airframe. 


Most of the degrees of freedom are 
reduced from the airframe gross stiffness 
matrix. Mass properties are concentrated 
at the remaining (retained) degrees of 
freedom. Equations (1) are the airframe 
equations of motion with the gross stiff- 
ness matrix. Equations (3) are the air- 
frame equations of motion, in terms of 
the reduced stiffness matrix. 
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Substructures Method 


A large saving in computer time can 
be realized by performing the matrix 
reduction process on several smaller sub- 
structure stiffness matrices instead of 
on the large stiffness matrix for the 
entire airframe. In one application, use 
of the substructures method reduced com- 
puter running time from about ten to two 
hours on an IBM 360/65 computer. 

The airframe is divided into several 
substructures, and all but mass and 
boundary degrees of freedom are reduced 
from the stiffness matrix of each sub- 
structure, The stiffness matrices of the 
substructures are then merged (super- 
imposed or added just as they are for 
individual members) to form a stiffness 
matrix for the entire airframe. Any 
degree of freedom on the boundaries may 
be reduced after merging the substructure 
matrices (Figure 3) . 

NAS TRAN 

New developments in finite element analy- 
sis have been occurring on a continuous 
basis. New programs and new structural 
elements, both dynamic and stress analy- 
sis capability, FORTRAN programming cap- 
ability by the engineer within the finite 
element program, and greater problem size 
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(9) 


capability have been developed^. NASTRAN 
(NASA Structural Analysis) 6 is a govern- 
ment developed, maintained, and continu- 
ally updated finite element program which 
has apparently provided a solution to the 
difficulties of developing and maintain- 
ing finite element programs by private 
contractors. NASTRAN is similar to other 
finite element computer programs except 
that it generally provides additional 
capability: 

• More types of structural elements 

• Common deck for stress and 
dynamic analysis 

w User programming capability 

• Transient vibration analysis, 
buckling, non-linear, and static 
capability 

• Unlimited size capability for 
mass and stiffness matrices. 


K} - _ M M hi 

The corresponding motions at the rotor 
hub are 

M = [ a RR“ a RA a aa a ar ] N ( 

The mobility matrices in the above equa- 
tions may be obtained analytically using 
computed modal properties (Equations (1) 
through (7)) or by applying unit vibra- 
tory loads to the airframe in a series of 
shake tests. 

This method was applied to predic- 
tion of cockpit vibration with a vertical 
cockpit absorber for the Model 347 heli- 
copter®. Analytical and flight test re- 
sults are compared in Figure 5 \ 

Coupled Rotor/Airframe Analysis 

Theory 


For a nominal fee, this program and 
manuals describing the program and its 
use are available. NASTRAN provides a 
standard for airframe dynamic analysis 
and relieves contractors of some of the 
problems of maintaining the most up-to- 
date methods for airframe structural 
analysis. 

Energy Methods for Structural Optimization 

One further development related to 
airframe dynamics is the Damped Forced 
Response Method 7 * 8. The airframe forced 
response is computed, and structural 
members with significant strain energy 
are identified. These members are 
changed to reduce vibration response for 
modes with frequencies above and below 
the rotor exciting frequency. This 
method is outlined in Figure 4. 

Vibration Control Devices 


Any vibratory motion of the rotor 
hub will change the rotor blade vibratory 
aerodynamic and inertia forces which are 
summed to obtain vibratory hub loads. 
Changes in hub loads will in turn cause 
changes in vibratory hub motions® 


Airframe Motion is assumed to be 
related to vibratory hub loads by a 
vnnbi 1 i l-v matrix for a particular fixcitina 


frequency: 


x ks 

< 


A kll A kl2 

* 

F ks 

x kc 


A k21 A k22_ 


F kc 


where 


-M 


F ks 

F kc 


(ID 


{ x k}={ X ks} sin knt + { X kc} cos knt 
{ F k} = { F ks} sin knt + ( F kJ cos knt 


Vibration control devices such as 
absorbers are often used to reduce vibra- 
tion in local areas of the airframe. 

The force output for an absorber may be 
computed by expressing the vibration as 
the sum of vibration due to rotor forces 
and the vibration due to the force output 
by the absorber. 




a aa 

a ar 


< 

[ X R 


a ra 

a rr 


f r 


The airframe mobility data are air- 
frame responses to unit vibratory hub 
loads? these data may be obtained analy- 
tically by using theoretical modal proper- 
ties (Equation (4) through (7)), or by 
conducting an airframe shake test. It is 
emphasized that these are airframe 
response characteristics for no blade 
mass attached to the airframe at the 
rotor hub. All blade inertia effects 
will be included in the rotor vibratory 
hub loads as modified by vibratory hub 
motion. 


The absorber force output required to null 
vibration at the absorber attachment 
point is 


In general, six sine and six cosine 
components of shaking forces and moments 
exist at each rotor hub? a tandem rotor 
helicopter would have a total of 24 
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components of vibratory forces. If only 
the rotor hub motions are considered, the 
relationship between hub motion and hub 
forces is: 


24x1 

24x24 

24x1 

X Hks 

35 

11 

F ks 

X Hkc 

i- J 

F kc 


( 12 ) 


The vibratory hub loads are assumed 
to be loads with no hub motion plus an 
increment of hub loads proportional to hub 
motion: 
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(13) 


The coefficients of the B matrix are 
obtained by making several computations of 
vibratory hub loads: 


• Components of vibratory hub loads 
are computed assuming no hub 
motion 

• Components of vibratory hub loads 
are computed assuming a small 
vibratory hub motion at the fre- 
quency for each degree of freedom 
of hub motion at each rotor 

• Changes in sine and cosine com- 
ponents of vibratory hub forces 
per unit vibratory hub motion in 
each rotor hub degree of freedom 
are then computed. 

The coupled rotor/airframe solution 
for compatible rotor hub motions and rotor 
hub loads is obtained by substituting 
Equation (13) in Equation (12) and solving 
for vibratory hub motions: 


x Hks 

, - 

[I]-[A] [B] 

pi 

[A], 

F kso 

x Hkc 




F kco 


Once a solution for Equation (14) is 
obtained, the total vibratory hub loads 
may be computed using Equation (13) and 
the vibration for the entire airframe may 
be computed using Equation (11) . 


Single Rotor Example 

Figure 6 shows a simple example of 
the coupled rotor airframe method applied 
to a single rotor helicopter vertical 
vibration analysis. Hub vertical vibra- 
tion response and the vertical vibratory 
hub loads are computed at a frequency of 
four times rotor speed (4/rev) . The air- 
frame is represented by its rigid body 
vertical mode and one flexible mode. 
Figure 6b shows airframe mobilities vs 
flexible mode natural frequency for 4/rev 
vertical hub forces. Hub vertical shak- 
ing forces vs hub vertical motion are 
shown in Figure 6c. The vibratory hub 
loads are seen to vary approximately 
linearly at least up to .005 inches of 
motion at the 4/rev frequency. Figures 
6d and 6e show compatible rotor hub 
vertical vibration amplitudes and rotor 
hub shaking forces vs flexible mode 
natural frequency . 

For this example, the rotor vibra- 
tory hub motions and forces both peak 
when the flexible mode natural frequency 
is just above the rotor hub force excit- 
ing frequency. This is not a general 
result, but depends upon the relation- 
ships between hub shaking forces and hub 
motions . 

Coupled Rotor/Airframe Analysis Computer 
PTogram (D-65) * 

Figure 7 shows the flow-diagram for 
the Boeing Vertol D-65 Coupled Rotor/ 
Airframe Analysis computer program. This 
program links three major computer 
programs 10 : 

Trim analysis program A-97 

Rotor vibratory hub loads analy- 
sis program D-88 

Airframe forced response analy- 
sis program D-96. 

Compatible fuselage motions and 
vibratory hub loads are obtained using 
this program with the method discussed 
above. In its current state, the D-65 
program computes three vibratory rotor 
forces and three vibratory rotor moments 
at each rotor for either single or tandem 
rotor helicopters. Response to trans- 
lational and rotational vibratory hub 
forces is computed for the airframe, but 
compatibility of hub forces and motions 
is satisfied for hub translational 
degrees of freedom only in the current 
version of the program. The program will 
be modified in the near future to provide 
compatibility for hub rotational degrees 
of freedom. 
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D.C., 1954. 


Analysis vs Test Results for the Model 347 
Helicopter 

The D-65 coupled rotor/airframe pro- 
gram was used to predict Model 347 flight 
vibration levels. Figure 8 shows the 
model used to predict airframe dynamic 
characteristics. Figure 9 compares pre- 
dicted vertical and lateral cockpit 
vibration levels vs vibration levels 
measured in flight. Vertical vibration 
levels are in reasonably good agreement 
at high airspeeds where vibration levels 
may become significant. Lateral vibra- 
tion levels are higher than predicted. 

Conclusions 


Methods have been developed indepen- 
dently for prediction of rotor vibratory 
hub loads and airframe dynamic character- 
istics. Methods are available for in- 
cluding effects of vibration control 
devices on airframe vibration and for 
optimizing the airframe structure. The 
substructure method is available for 
minimizing computer running time in 
analysis of airframe structures, and 
NASTRAN now provides a common finite 
element structural analysis program avail- 
able to all aerospace contractors. Rotor 
hub vibratory motions can modify rotor 
hub vibratory forces acting on the air- 
frame. A linear coupled rotor/airframe 
analysis method provides an approach for 

1 4 * — ‘*-11 1 . . ^ ^ J 1 , 1 

U-C W L4. XU- 14. L&Jy/ 1UV/ WAW14W m 

shaking forces. This method should be 
studied further to determine its 
validity. A method of this type should 
be considered in applications of NASTRAN 
for prediction of helicopter vibration; 
the user programming feature in NASTRAN 
should permit a coupled rotor/airframe 
solution of this type within NASTRAN. 

Figure 10 shows a scheme for solving 
for rotor trim, rotor forces with no hub 
motion, and the rotor impedance matrix 
using a rotor analysis program. NASTRAN 
would be programmed to use these mobili- 
ties and the rotor analysis results to 
solve for compatible rotor/airframe loads 
and motions. The NASTRAN airframe analy- 
sis could include airframe installed 
vibration control devices either in the 
initial airframe analysis or in the 
coupled rotor/airframe solution. Finally, 
results of these analyses could be used 
to determine optimum changes to the air- 
frame structural' elements for minimizing 
airframe vibration. 
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Figure 1. Rotor Blade Analytical Model 



Figure 2. Uncoupled Airframe Dynamic 
Analysis 



Figure 2 . Continued 
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OPTION 1 


• AIRFRAME INPUT 
DATA, SUBSTRUCTURE i 

i 

• GENERATE SUBSTRUCT. 
i GROSS STIFFNESS 
MATRIX 

• REDUCE ALL BUT MASS 
AND BOUNDARY DEGREES 
OF FREEDOM 


• MERGE SUBSTRUCTURE REDUCED 
STIFFNESS MATRICES 


• REDUCE NON-MASS BOUNDARY 
DEGREES OF FREEDOM 


L 

1 

• COMPUTE AI 
MODES AND 

RFRAME NATURAL 
FREQUENCIES 


• ROTOR FORCES 


• COMPUTE AIRFRAME 
FORCED RESPONSE 



B 




UJJ 

— BOUNDARY DEGREES 
OF FREEDOM 




Figure 3, Substructure Method for Gener- 
ating Airframe Reduced Stiff- Figure 4. Damped Forced Response Method 

ness Matrix for Airframe Optimization 


87 






o 

+ 1 .4 h 


§ 3 

H ' J 


2 *5 1 / , 

« - 1 7 - 


« 

u 


MODEL 347 COCKPIT VERTICAL 
(NO ABSORBERS) 

Z H4s 

SB 

A 11 a 12 


f Z4s 

FLIGHT DATA p 

_ Z H4c 


A 21 A 22 


F Z4c 


a 22 = A 11 


A 21 A 12 


OBJECTIVE 

ANALYSIS 



40 80 120 160 

AIRSPEED - KNOTS 
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(b) Airframe -Hub Mobilities 


Figure 5, Predicted Vs . Measured 

Cockpit vibration Reduction 
with a Vertical Cockpit 
Absorber 
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(a) Single Rotor Helicopter Vertical 
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Vibration 


(c) Hub Forces Vs. Hub Motion 


Figure 6. Coupled Rotor /Air frame Analysis 
for a Single Rotor Helicopter 
Vertical Vibration 


Figure 6. Continued 
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(d) Vibratory Hub Motion 
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(e) Vibratory Hub Force 
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3. G-Levels 


Figure 6. Continued 


Figure 7. D— 65 Coupled Rotor /Air frame 
Program Flow Diagram 
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Figure 8. Model 347 Airframe Dynamic Model 
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Figure 9. Model 347 Flight Data Vs. D-65 
Coupled Rotor/Airframe Analysis 
Results 


• TRIM ANALYSIS 


• VIBRATORY ROTOR LOADS , 
NO HUB MOTION {f r0 } 

• VIBRATORY ROTOR LOADS 
WITH UNIT VIBRATORY 
HUB MOTIONS 


• ROTOR IMPEDANCE 
MATRIX, B 


• AIRFRAME SUBSTRUCTURE ANALYSIS 
(INCLUDE MODELS OF VIBRATION 
CONTROL DEVICES) 

•MERGE SUBSTRUCTURE STIFFNESS 
MATRICES 

•COMPUTE- AIRFRAME MODES , FREQUENCIES , 
AND GENERALIZED MASSES WITH NO BLADE 
MASS AT ROTOR HUBS 


•COMPUTE AIRFRAME RESPONSE TO UNIT 
VIBRATORY HUB LOADS: 
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• COMPUTE 
MOTIONS 

COMPATIBLE VIBRATORY HUB 
AND FORCES 

, ~~sr 

II 

i] “ [ a rr] [ B l] [ a rr] { f ro) 

• COMPUTE 

TOTAL HUB FORCES 

N ■ 

M + HM 

• COMPUTE 
DEGREES 

MOTIONS AT OTHER AIRFRAME 
OF FREEDOM 

M - 

[ a ar] ( f r] 
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• IDENTIFY STRUCTURAL CHANGES 
TO MINIMIZE AIRFRAME 
VIBRATION USING STRAIN 
ENERGY METHODS 


Figure 10. Coupled Rotor/Airf rame/NASTRAN 
Analysis 
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